Command line: [exonerate -m p2g --showtargetgff -q ./gpx3human.fa -t ./gpx3human/subseq/gi|432098791|gb|KB109247.1|_subseq.fa] Hostname: [luke] C4 Alignment: ------------ Query: SPP00000021_2.0 # Protein # Glutathione peroxidase 3 (GPx3) # Human (Homo sapiens) Target: gi|432098791|gb|KB109247.1|:subseq(1939698,100182) Myotis davidii unplaced genomic scaffold scaffold18, whole genome shotgun sequence Model: protein2genome:local Raw score: 916 Query range: 28 -> 226 Target range: 48237 -> 51006 29 : LysMetAspCysHisGlyGlyIleSerGlyThrIleTyrGluTyrGlyAlaLeuThrIle : 48 :!!! !|||||||||!.!|||:!:|||||||||||||||||||||||||||||||||:!! GlnThrAspCysHisAlaGlyValSerGlyThrIleTyrGluTyrGlyAlaLeuThrLeu 48238 : CAGACGGACTGCCATGCTGGTGTGAGCGGCACCATCTATGAGTATGGAGCCCTCACCCTC : 48295 49 : AspGlyGluGluTyrIleProPheLysGlnTyrAlaGlyLysTyrValLeuPheValAsn : 68 :!!||||||||||||||||||||||||||||||||||||||||||||||||||||||||| AsnGlyGluGluTyrIleProPheLysGlnTyrAlaGlyLysTyrValLeuPheValAsn 48296 : AACGGGGAGGAGTACATCCCCTTCAAGCAGTATGCCGGCAAATACGTCCTCTTTGTCAAC : 48355 69 : ValAlaSerTyrUnkGlyLeuThrGlyGlnTyrIle{G} >>>> Target Intron : 81 |||||||||||| |||||||||.!!||||||:!!{|} 1110 bp ValAlaSerTyr***GlyLeuThrSerGlnTyrVal{G}++ 48356 : GTGGCCAGCTACTGAGGCCTGACGAGCCAGTACGTT{G}gt................... : 48397 82 : 1 >>>> {lu}LeuAsnAlaLeuGlnGluGluLeuAlaProPheGlyLeuValIleLeu : 97 {||}|||||||||||||||||||||||| !||||||||||||||||||||| ++{lu}LeuAsnAlaLeuGlnGluGluLeuGlnProPheGlyLeuValIleLeu 48398 : ......ag{AA}CTGAATGCACTACAGGAAGAGCTGCAACCATTTGGTCTGGTCATTCTG : 49552 98 : GlyPheProCysAsnGlnPheGlyLysGlnGluProGlyGluAsnSerGluIleLeuPro : 117 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| GlyPheProCysAsnGlnPheGlyLysGlnGluProGlyGluAsnSerGluIleLeuPro 49553 : GGCTTCCCCTGCAACCAATTTGGAAAACAGGAGCCAGGAGAGAACTCGGAGATTCTACCT : 49612 118 : ThrLeu{Ly} >>>> Target Intron 2 >>>> {s}TyrValArgProGlyGly : 126 :!!|||{||} 445 bp {|}:!!||||||||||||||| SerLeu{Ly}++ ++{s}HisValArgProGlyGly 49613 : TCCCTC{AA}gt.........................ag{G}CATGTTCGACCAGGTGGG : 50084 127 : GlyPheValProAsnPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGluGln : 146 |||!:!|||||||||||||||||||||||||||||||||||||||||||||!!.|||||| GlyTyrValProAsnPheGlnLeuPheGluLysGlyAspValAsnGlyGluAsnGluGln 50085 : GGCTACGTCCCCAATTTCCAACTCTTTGAGAAAGGAGATGTGAACGGGGAGAACGAGCAG : 50144 147 : LysPheTyrThrPheLeuLys >>>> Target Intron 3 >>>> AsnSerCysP : 157 ||| !!||||||||||||||| 620 bp |||||||||| LysValTyrThrPheLeuLys++ ++AsnSerCysP 50145 : AAGGTCTACACGTTCCTGAAGgt.........................agAACTCCTGCC : 50797 158 : roProThrSerGluLeuLeuGlyThrSerAspArgLeuPheTrpGluProMetLysValH : 177 ||||||||!!!||||||||||||:!! !! ! ||||||||||||||||||||||||:!!| roProThrThrGluLeuLeuGlySerProLysArgLeuPheTrpGluProMetLysIleH 50798 : CTCCCACCACGGAGCTCCTGGGCTCGCCTAAACGCCTCTTCTGGGAACCCATGAAGATCC : 50857 178 : isAspIleArgTrpAsnPheGluLysPheLeuValGlyProAspGlyIleProIleMetA : 197 |||||||||||||||||||||||||||||||||||||||||||||||! !|||||||||| isAspIleArgTrpAsnPheGluLysPheLeuValGlyProAspGlyLysProIleMetA 50858 : ATGACATCCGCTGGAACTTTGAGAAGTTCCTGGTGGGGCCGGACGGTAAACCCATCATGC : 50917 198 : rgTrpHisHisArgThrThrValSerAsnValLysMetAspIleLeuSerTyrMetArgA : 217 |||||:!!|||||||||||||||||||||||||||||||||||||||:!!|||:!!|||. rgTrpTyrHisArgThrThrValSerAsnValLysMetAspIleLeuAlaTyrLeuArgG 50918 : GCTGGTACCACCGGACCACGGTCAGCAACGTCAAGATGGACATCCTGGCCTACTTGAGGG : 50977 218 : rgGlnAlaAlaLeuGlyValLysArgLys : 226 .!||||||||||||!.!|||:!!..!||| luGlnAlaAlaLeuAlaValGluGluLys 50978 : AGCAGGCAGCCCTGGCGGTCGAAGAGAAG : 51006 vulgar: SPP00000021_2.0 28 226 . gi|432098791|gb|KB109247.1|:subseq(1939698,100182) 48237 51006 + 916 M 52 156 S 0 1 5 0 2 I 0 1106 3 0 2 S 1 2 M 38 114 S 0 2 5 0 2 I 0 441 3 0 2 S 1 1 M 33 99 5 0 2 I 0 616 3 0 2 M 73 219 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2013-11-27 ##type DNA # # # seqname source feature start end score strand frame attributes # gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local gene 48238 51006 916 + . gene_id 1 ; sequence SPP00000021_2.0 ; gene_orientation + gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local cds 48238 48394 . + . gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local exon 48238 48394 . + . insertions 0 ; deletions 0 gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local splice5 48395 48396 . + . intron_id 1 ; splice_site "GT" gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local intron 48395 49504 . + . intron_id 1 gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local splice3 49503 49504 . + . intron_id 0 ; splice_site "AG" gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local cds 49505 49622 . + . gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local exon 49505 49622 . + . insertions 0 ; deletions 0 gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local splice5 49623 49624 . + . intron_id 2 ; splice_site "GT" gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local intron 49623 50067 . + . intron_id 2 gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local splice3 50066 50067 . + . intron_id 1 ; splice_site "AG" gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local cds 50068 50167 . + . gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local exon 50068 50167 . + . insertions 0 ; deletions 0 gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local splice5 50168 50169 . + . intron_id 3 ; splice_site "GT" gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local intron 50168 50787 . + . intron_id 3 gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local splice3 50786 50787 . + . intron_id 2 ; splice_site "AG" gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local cds 50788 51006 . + . gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local exon 50788 51006 . + . insertions 0 ; deletions 0 gi|432098791|gb|KB109247.1|:subseq(1939698,100182) exonerate:protein2genome:local similarity 48238 51006 916 + . alignment_id 1 ; Query SPP00000021_2.0 ; Align 48238 29 156 ; Align 49507 82 114 ; Align 50069 121 99 ; Align 50788 154 219 # --- END OF GFF DUMP --- # -- completed exonerate analysis